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ABSTRACT 



We study four scenarios for the SCP 06F6 transient event that was announced 
recently. Some of these were previously briefly discussed as plausible models for 
SCP 06F6, in particular with the claimed detection of a. z = 0.143 cosmological 
redshift of a Swan spectrum of a carbon rich envelope. We adopt this value of z 
for extragalactic scenarios. We cannot rule out any of these models, but can rank 
them from most to least preferred. Our favorite model is a tidal disruption of a 
CO white dwarf (WD) by an intermediate- mass black hole (IMBH). To account 
for the properties of the SCP 06F6 event, we have to assume the presence of a 
strong disk wind that was not included in previous numerical simulations. If the 
IMBH is the central BH of a galaxy, this explains the non detection of a bright 
galaxy in the direction of SCP 06F6. Our second favorite scenario is a type la-like 
SN that exploded inside the dense wind of a carbon star. The carbon star is the 
donor star of the exploded WD. Our third favorite model is a Galactic source of 
an asteroid that collided with a WD. Such a scenario was discussed in the past 
as the source of dusty disks around WDs, but no predictions exist regarding the 
appearance of such an event. Our least favorite model is of a core collapse SN. 
The only way we can account for the properties of SCP 06F6 with a core collapse 
SN is if we assume the occurrence of a rare type of binary interaction. 



1. INTRODUCTION 

SCP 06F6 is an optical transient discovered by Barbary et al. (2008). It raised to an 
optical peak flux of Fpeak — 2.5 x lO"^"^ erg s~^ cm~^ in ~ 100 days, and then dechned in 
another ~ 100 days. Barbary et al. (2008) ruled out a microlensing event, and noted that if 
SCP 06F6 is associated with the galaxy cluster CL 1432.5+3332.8 at a redshift of z = 1.112, 
then its peak luminosity is like that of the most bright SNe observed to date. In this case. 
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this object projected distance from the center of the cluster would be 290 kpc, but (1) there 
is no obvious host galaxy, making the association questionable (a faint host is possible at a 
projected distance of 12kpc for this redshift; Barbary et al. 2008). (2) The spectrum of this 
object is not like any other SN. 

Gaensicke et al. (2008a) presented a fit of the spectrum of SCP 06F6 with a spectrum 
of a carbon-rich star at a redshift oi z — 0.143. They also argue that the hghtcurve looks 
somewhat like that of a supernova type II. Based on that they suggest SCP 06F6 represents 

a new class of core collapse SNe. However, this explanation shares the weak points (1) and 
(2) above, and in addition: (3) The X-ray luminosity reported by Gaensicke et al. (2008a) 
is two orders of magnitude above that of a typical SN. (4) The spectral fit is not perfect. 

We discuss alternative possible scenarios for the SGP 06F6 transient. In section 2 we 
describe the critical observations that set constraints on any model to explain the transient. 
In section 3 we discuss a Galactic WD-asteroid merger model. Barbary et al. (2008) consid- 
ered as well the possibility that the progenitor was a cool WD in our galaxy. In that case the 
WD progenitor should be quite cold, and its distance should he D > 1.2 kpc, assuming the 
WD is not colder than 3000 K. In section 4 we discuss three possible extragalactic models. 
In addition to the core collapse SN model (Gaensicke et al. 2008a) we discuss a type la-hke 
SN scenario, where the exploding WD is enshrouded in a carbon-rich nebula formed by the 
wind from the WD's giant companion. We also discuss a possible disruption of a WD passing 
close to an intermediate- mass black hole (IMBH; Rosswog et al. 2008b, c). We compare the 
four models with observations and summarize in section 5. 



2. CRITICAL OBSERVATIONS 

A model aimed at explaining SGP 06F6 should be able to simultaneously explain its 
different observational features. The crucial constraints are as follows (see also Table 1). 

(1) The Hghtcurve. The lightcurve of SGP 06F6 is symmetric in time and akin to a bell 
curve. The duration of the transient is ~ 200 days, and only near the peak it resembles the 
lightcurve of a typical SN. At later times SNe decline at a much lower rate than the late 
evolution of SGP 06F6. 

(2) The spectrum. Barbary et al. (2008) checked that the SGP 06F6 spectrum does 
not resemble any known SN spectra. As the closest, but not very convincing, matches they 
suggested broad absorption hue quasars or white dwarfs (WDs) with carbon hues (DQ WDs). 
Indeed, the shapes and spacing of the absorption features in SGP 06F6 bear resemblance to 
the G2 Swan bands, present in DQ WDs and carbon stars. 
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Determination of temperature from the spectrum partially depends on the object's 
galactic or extragalactic origin. Barbary et al. (2008) estimated Tes > 6500 K for a galactic 
object from the slope of the red continuum. Gaensicke et al. (2008a) hold that the presence 
of Swan bands implies a temperature of 5000 - 6000 K. However, Swan bands are observed 
at higher temperatures, e.g. in DQ WDs at up to ~ 10,000 K (Dufour et al. 2005). They 
generally weaken with increasing temperature, but in a carbon-dominated atmosphere they 
would be extremely strong even at 10,000 K (Dufour et al. 2008). 

Barbary et al. (2008) find some evidence for spectral evolution in the IR colors and the 
three optical spectra, but it is only mild. The photosphcric temperature remains apparently 
almost constant over the 36 day period spanned by the three optical spectra. The two main 
model constraints coming from the spectrum arc a relatively cool photosphere (~ 6000K) 
and the presence of very wide bands, most likely C2 Swan bands. 

(3) X-ray emission. The X-ray flux was measured at one point in decline, about 80 days 
after maximum, when the optical flux was at ~ 20 — 30% of its peak value. The X-ray flux 
was ~ 4 times the peak optical flux, i.e., ~ 10 — 20 times the optical flux at the time. It is 
quite possible that most of the emission from SCP 06F6 came in the X-ray band. 

The presence of both X-rays and a photosphere at T < 6500 K can be used to constrain 
the geometry of the event. Most likely, the X-rays result from post-shocked gas that was 
ejected by the progenitor at velocities of several x 1000 km s~^. To preserve the photosphere 
at low temperatures, though, the X-ray emitting gas has to be ejected over a small solid angle. 
The most hkely geometry is that of jets launched in the polar directions, with the observer 
not along the polar axis; we rather face the disk sides. Such a geometry indicates large 
specific angular momentum, which in turn promotes models that involve binary systems, 
e.g., merger, tidal disruption, or binary progenitor model. 

3. GALACTIC ORIGIN: A WD-ASTEROID COLLISION 

Although the galactic origin model is not our favorite one, we discuss it below for three 
reasons. (1) We cannot completely rule it out, and for that it serves us also in discussing 
the other models. (2) It might be applicable to other transient events in the future. (3) 
It draws some parallel ingredients with the merger model of V838 Mon. In addition, the 
process of WD-asteroid collision and the formation of circum-WD material was discussed 
before in other contexts, but no predictions exist regarding the appearance of such an event. 

Some WDs are thought to possess dusty disks because of their observed infrared excess 
(e.g. Zuckerman & Becklin 1987; Graham et al. 1990; Becklin et al. 2005; Kilic et al. 
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2006; von Hippel et al. 2007). The evidence was strengthened by the discovery that these 
WDs also have metal-rich photospheres, presumably due to accretion of metal-rich material 
(Koester et al. 1997; Zuckerman & Reid 1998; Jura et al. 2007; Zuckerman et al. 2007). 
Even more recently, some hotter WDs have been found which exhibit double-peaked Ca II 
emission lines while at the same time lacking Balmer emission (Gaensicke et al 2006, 2007, 
2008b). The double-peaked Ca II emission profiles allow to conclude that the metal-rich 
gas around these WDs has indeed disk-like geometry. As discussed by Debes & Sigurdsson 
(2002) and Jura (2003), such disks can result from a tidal disruption of an asteroid or a 
comet colliding with the WD. In this scenario, asteroids or comets are thought to be sent on 
a collision course with the WD by some perturbation in the asteroid belt or cometary cloud 
surrounding the WD. 



3.1. Model ingredients 

Consider the well studied transient event V838 Mon (e.g., Wisniewski et al. 2008; 
Kaminski 2008; Sparks et al. 2008; and papers in Corradi & Munari 2007). The most 
popular model for the eruption of V838 Mon is a merger event, termed mergeburst, of two 
main sequence (or pre-MS) stars (Soker & Tylenda 2003, 2006; lyienda 2005; Tylenda & 
Soker 2006). The basic ingredients of the merger model are as follows: (1) Two objects 
collide, and the less dense object is completely destroyed (Shara 2002). (2) The collision 
sometimes creates Shockwaves that might lead to X-ray emission (Shara 2002). (3) Some 
fraction of the mass of the less dense object is accreted within a short time by the denser 
object. This is the source of the large energy in the flash of the event. (4) The other fraction 
of the mass of the destroyed objects, and possibly a small amount of mass from the massive 
object, are ejected to large distances, and form an extended envelope. Some fraction of the 
extended envelope might escape the system. (5) The release of gravitational energy continues 
as the inner parts of the extended envelope contract and are accreted (Tylenda 2005). (6) 
Because the release of the gravitational energy is due to a contracting envelope, the energy 
release will continue as long as there is an envelope. For that, in the bright phase we will 
not sec the hot inner parts (assuming spherical symmetry; see below). The merger will fade 
as a relatively cool object. This is in contrast to novae, where the envelope disperses while 
the inner surface of the hot WD is still hot: novae get bluer as they fade. 

V838 Mon transient event is different in its spectrum, lightcurve, and timescale from 
SCP 06F6. However, we will examine what is required from a merging Galactic object to 
explain an event like SCP 06F6. Our proposed model is based on the same physical processes 
as in the model for V838 Mon, replacing two MS stars with a large asteroid hitting a WD. 
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3.2. Energy considerations 

Let the distance to SCP 06F6 he D = lODio kpc. Its peak luminosity is Lp = 
O.OSDIqLq, and the total radiated energy in the optical bands is (assuming the source is 
in our galaxy) 

^rad ^ 3 X 10-^^D?o erg (1) 

For a black body temperature of 6500 K the photospheric radius of the object at peak 
luminosity is 

i?ph ~ 0.2D,oRq. (2) 

Including the X-ray emission, the total radiated energy might be much larger, by a factor of 
~ 5 (Gaensicke et al. 2008a). 

If the energy is released by accreting an amount of mass m^, then the energy released 
to radiation and to the inflated envelope is E = GMyf-£,ma/2R-ysiD. Let a fraction rj of this 
energy be radiated in the optical bands (by the photosphere). From their study of the stellar 
collision model to the V838 Mon transient event, Soker & Tylenda (2003) estimate that 
7] ~ 0.025—0.1. The rest of the gravitational energy will be released later, after the envelope is 
dispersed. The accreted mass based on the observations is then ma ~ 5 x IO^^D^qT/"^ g, where 
we took for the WD mass and radius Mwd = 0.6Mq and Rwd — 8, 000 km, respectively. If 
half the mass is accreted, and half resides in the envelope, then the total mass of the object, 
an asteroid or a comet, is 

m, ^ 10^' Dl[^y' g. (3) 
This is an asteroid with a diameter of ~ 1000 km. 



3.3. Extended envelope structure 

The observations indicate a black body temperature of ~ 6500 K (Barbary et al. 2008). 
We note that at this temperature the opacity is /t ~ 0.1 cm^ g^^, and increases with tem- 
perature. The dependence of opacity on temperature in the ~ 4, 000 — 10, 000 K region 
is K ~ 10(T/10'^ K)^° cm^ g"^ (Alexander & Ferguson 1994; Marigo 2002; Ferguson et al. 
2005). The average density in the envelope is p ~ 3 x 10~^D^q^(?7/0.1)"^ g cm~^. The optical 
depth is T ~ nRphP — 50(77/0.1)"^. So an optically thick envelope is formed. 

This actually sets the typical radius of the envelope for such kind of objects. We demand 
an optical depth of r ~ 1 at the photosphere, where the density is much lower than the 
average density. Therefore, the average optical depth of the envelope will be much higher. 
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and we take Tg ~ 10. For an envelope mass of 5 x 10^^(?7/0.1) ^ g and k — 0.1 cm^ g ^, the 
condition Tg ~ reiilphP = 10 gives i?ph = 0.5(^/0.1 cm^ g-i)V2^-V2(^^/io)-V27^Q. 

The effective temperature goes as Tph oc R^^^'^L]Jp , where Lap is a possible alternative 
luminosity. Using this expression in the dependence of opacity on temperature, and then in 
the condition r = Tg, gives 

The energy released is fixed by the asteroid that collided with the WD. So the typical 
luminosity (about the peak luminosity) depends on the time duration of burst th as Lap oc t^^. 
Hence the radius depends on the duration of the event as i?ph oc t^°'^^. Namely, if the event 
was longer (shorter) by a factor of 10, the radius would be smaller (larger) by a factor of 2.3. 
The dependence on the exact distribution of energy in the event, as expressed by 77, is also 
very weak. Overall, we expect that for the WD and asteroid used here, the photospheric 
radius would be i?ph — 0.1 — li?©, with a very weak dependence on the unknown variables. 

The same arguments show that the steep increase of the opacity with temperature means 
that there will be no strong temperature evolution with time. Again, the reason is that as 
the envelope thins, the photosphere gets deeper to hotter layers. But there the opacity 
increases, so the difference in the physical depth into the envelope is not large. 

The Kelvin- Helmholtz time of the envelope is 

tKH- -:.100i^, (^j days (5) 

Overall, we can account for the optical lightcurve by emission from a semi-hydrostatic 
envelope that is dispersed in a Kelvin-Helmholtz timescale. The envelope is not in full 
hydrostatic equilibrium. If it deviates from a static configuration with a velocity of ~ 0.1 
times the photospheric sound speed, then the dispersion occurs over about one month. Most 
of the envelope mass will be accreted onto the WD. But now, as we can see deeper into 
the WD, and because of the angular momentum of the colliding asteroid, the system will 
resemble a dwarf nova in an outburst. 

At this stage we expect strong X-ray emission as well. This is a possible explanation for 
the strong X-ray emission from SCP 06F6 in the WD-asteroid merger model. If indeed most 
of the emission comes in the X-ray band, our explanation for that is the emission from near 
the surface of the WD. As the thermal energy content of gas near the surface of the WD is 
larger, the Kelvin-Helmholtz time will be somewhat longer for the X-ray emission and the 
optical phase will be followed by an X-ray luminous phase. 
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The study of the X-ray emission requires 3D hydrodynamical simulation of the merger 
event. However, we do note that accreting WDs can have a large fraction of their emission in 
the X-rays. For example, RU Peg has ~ 1/3 of its radiation coming in X-rays in quiescence, 
when its accretion rate is ~ 3 x 10^^ g s^^ (Silber et al. 1994). The accretion rate during 
the decline of the WD-asteroid model is similar. Some other WDs accreting in quiescence 
have their X-ray emission about equal to that in the optical (van Teeseling et al. 1996). The 
condition for that is that we see the WD surface. This can occur in our merger model if 
after the optical peak the envelope becomes clumpy, or because of the angular momentum 
of the asteroid, a very oblate atmosphere is formed, which later flattens into a disk. 

We note that equation 1^ ruled out a WD- WD merger. A WD- WD merger was studied 
in the past in relation to the formation of R Coronae Borealis (RGB) stars (e.g., Iben et al. 
1996), and in principle can be a source of a transient event. However, for the luminosity 
and radius inferred for SCP 06F6, the Kelvin-Helmholtz time is much too long. Also, a 
non-exploding WD- WD merger cannot be an extragalactic source at z = 0.143, because the 
required luminosity would be much larger than the Eddington luminosity of such a system. 

3.4. Spectrum 

In the WD-asteroid model the spectral features seen in SCP 06F6 are caused by the C2 
Swan bands (or molecules containing also hydrogen) at 2; = 0. The observed optical spectrum 
is not easy to explain in this the bands observed in SCP 06F6 do not exhibit the 

typical placement and shapes observed in Swan bands in most objects. However, distortions 
and displacements of the C2 Swan bands are observed among the DQ WDs, especially in 
the peculiar DQp type (Hall & Maxwell 2008). Their exact cause is still disputed, the main 
candidates being magnetic fields or molecules related to, but different from C2, such as C2H 
(see Schmidt et al. 1995, 1999; Hall & Maxwell 2008). Figure [T] compares the spectrum 
of SCP 06F6 taken by KECK on 28.05.2006 (Barbary et al. 2008) with spectra of various 
objects exhibiting Swan bands. With the exception of comet Halley (de Freitas Pacheco et 
al. 1988) and LP 790-29 (Schmidt et al. 1995) other spectra used for comparison are from 
the Sloan Digital Sky Survey (SDSS, York et al. 2000; Adelman-McCarthy et al. 2008). 

We note that the positions of the redder spectral features of SCP 06F6 coincide rea- 
sonably well with the positions of normal Swan bands (bandheads of the latter are marked 
by vertical dotted lines in Fig. [1]), even though the shapes tend to be more symmetric and 
rounded in the transient than in normal DQ WDs or carbon stars. The general shape of the 
bands in SCP 06F0 is more similar to the DQp stars than to the other objects featuring C2. 
To the blue, the broad maximum at 4500-4700A in SCP 06F6 is clearly discrepant with the 
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Fig. 1. — Comparison of the optical spectrum of SCP 06F6 (taken on 28.05.2006 with KECK, 
Barbary et al. 2008) to various objects exhibiting Swan bands. The spectra, normahzed to 
1 at 6500A and shifted vertically, are labeled with the categories and names of the objects. 
The top spectrum is of the same carbon star as found by Gaensicke et al. (2008a) to be, 
when redshifted to z=0.143, the best match for SCP 06F6. To assess their fit, this spectrum 
should be compared to the second spectrum from the top, which is the spectrum of SCP 06F6 
in its assumed restframe, if the transient is indeed located at z=0.143. Dotted vertical lines 
mark band-heads of the C2 Swan bands. Two short dashed vertical lines connect the features 
of the two DQp WDs with an otherwise unaccounted for maximum at 45OO-47OOA in the 
SCP 06F6 spectrum. Note that the comet spectrum (bottom curve) features Swan bands in 
emission. 
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normal Swan bands and also with most DQp WDs. It is only matched by a corresponding 
maximum in the spectra of peculiar magnetic DQp WDs LP 790-29 and SDSS J1113+0146 
(see the two short dashed vertical lines in Fig. [1]). 

A potential difficulty to the C2 Swan band explanation of the SCP 06F6 spectrum is that 
normally the central bands are the strongest (usually the Af = band with head at 5165A) 
and the bands to the sides weaken gradually. That is not the case for SCP 06F6, which has 
an alternating weak-strong band intensity pattern. In fact, this issue is not specific to the 
galactic model. It is then fortunate that we can point to a magnetic DQp WD with Swan 
bands and a similar band intensity pattern - namely G 99-37 (SDSS 0548-1-0010) which has 
B ~ lOMG (Angel 1977). In this case the alternating pattern results from a presence of a 
strong CH band at 4300A in addition to the C2 bands (e.g., Dufour et al. 2005). 

There is no one ideal fit, but the best case seems to be LP 790-29, a strongly magnetized 
WD with B ~ 50MG and T^s = 7800K (Liebert et al. 1978; Bues 1999). The lack of a perfect 
match is not surprising as the temperature of LP 790-29 is higher than that of the transient 
photosphere. Also, the magnetic field in SCP 06F6 may be somewhat stronger, since its blue 
maximum is more blueshifted than in LP 790-29, but less than in SDSS J1113-t-0146 (for 
which there is no B determination). 

Therefore, the variations in the C2 bands observed among DQp WDs could quite possibly 
account for the positioning of the SCP 06F6 features. The requirement is that the magnetic 
field in the photosphere, at a radius of i?ph — 20Z)io-Rwd, where -Rwd is a typical radius 
of a WD, be strong enough to modify the spectrum, i.e. ~ lOOMG. If the magnetic field 
originates at the WD surface, then even assuming magnetic field decreasing like 1/R^, the 
field at the WD surface would have to be 40,000D^qMG. This is a tall order, since the 
highest magnetic fields observed in WDs are ~ lOOOMG (Wickramasinghe & Ferrario 2000; 
Vanlandingham et al. 2005; Kawka et al. 2007). Still, it can be fulfilled if the object is 
relatively close, Diq < 0.16, not much more distant than the lower limit of 1.2kpc estimated 
by Barbary et al. (2008). Magnetic field amplification in the accretion disk may be another 
solution to account for the required high B. 

Alternatively, some misalignments (most notably the peak around 4500-4700A ) can 
be understood if strong emission is present on top of strong absorption bands. Swan bands 
emissions are well known in comets (e.g., Krishna Swamy 1997), they are also observed 
in R CrB stars in decline phase (Rao & Lambert 1993). Spectrum of comet Halley (from 
de Freitas Pacheco et al. 1988) is shown as an example in Fig. [H The cometary Swan 
bands emissions originate from resonant fluorescence (Stawikowski & Swings 1960) and their 
relative intensities vary between comets (e.g.. Fink & Hicks 1996) and within the same comet 
both spatially and temporally (O'Dell et al. 1988). The hypothetical Swan band emissions 



in SCP 06F6 need not be the same as in comets, but this example illustrates that some of 
the absorption features could be strongly affected (displaced or even reversed) while others 
still appear mostly unchanged. 

Figure [T] presents also a comparison with the spectral fit of Gaensicke et al. (2008a) 
who obtain z = 0.143 for SCP 06F6. The top curve displays their best matching carbon 
star spectrum, SDSS J0018-1101. Since all the comparison spectra in Fig. [T] are displayed 
at z=0, assessing the fit of Gaensicke et al. (2008a) requires also the transient spectrum 
to be transferred to its restframe. The second curve from the top represents the restframe 
spectrum of SCP 06F6, assuming the transient is indeed located at z=0.143. It can be noted 
that their fit is appealing but also not ideal, with each band having different deviations from 
their best case C star. 



/^From their spectral fit, Gaensicke et al. (2008a) determine a redshift of z = 0.143 for 
SCP 06F6 and, consequently, favor an extragalactic explanation. This redshift estimate being 
quite compelling, we adopt it for our exploration of the extragalactic scenarios. This sets 
immediate constraints for the energetics of the event and the size of the observed photosphere. 

The total peak observed flux derived from the ^775 band magnitude and from the spectral 
energy distribution (IR/visible/UV) is 2.5 x 10~^^ erg s~^ cm~^ (Barbary et al. 2008). At 
a redshift of z = 0.143 this translates to a peak luminosity of Lp ~ 1.3 x 10*^^ erg s~^. If 
the X-ray emission was in fact four times stronger, as it was later at the time of the X- 
ray observation (Gaensicke et al. 2008a), then Lp^x — 5 x 10^^ erg s~^. For a black body 
temperature T^g and spherical geometry, the radius of the photosphere at the peak luminosity 
is 



The lightcurve allows also a simple estimate of the total energy radiated through the duration 
of the event, tevent, as E^ad - -^^p4vent/2 ~ 10^^(Lp/10'^2 erg s"^) erg. 



Both Barbary et al. (2008) and Gaensicke et al. (2008a) consider SN origins for SCP 06F6. 
While if at z = 1.112 (as in Barbary et al.) the SN would have to be among the most lu- 
minous observed just to account for the optical flux, the z = 0.143 value does not require 
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4.1. Supernova Models 
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the conjectured SN to be over luminous. Typical Lp of a supernova is 10^^-10^^ erg s~^, with 
core collapse SNe being usually less bright but exhibiting larger brightness diversity than 
thermonuclear ones (e.g., Contardo et al. 2000; Richardson et al. 2002; Pastorello et al. 
2005). Typical radiated energy is 10^^ erg. Therefore the observed peak luminosity and en- 
ergy budget of SCP 06F6 are consistent with the supernova hypothesis, be it thermonuclear 
(like in SNe la) or core collapse. SNe around maximum light exhibit spectral evolution on a 
timcscale of days to weeks, changing the appearance of spectral lines and becoming gradually 
redder (e.g.. Wells et al. 1994; Filippenko 1997). Only slight spectral evolution was found 
by Barbary et al. (2008) for SCP 06F6, but given the scarce spectral data for the transient, 
there is no real inconsistency in this regard. 

SNe arc also known to be X-ray sources, but their X-ray luminosities arc 100 times lower 
than observed in this transient. However, we note that the kinetic energy of SNe ejecta is 
typically ~ 10^^ erg for both la and core collapse SNe (e.g. Wheeler et al. 1995) which is 
~ 100 times larger than that observed as radiation. Even higher kinetic energies are observed 
in the so called hypernovae, which are defined by their kinetic energy exceeding 10^^ erg 
(Paczyhski 1998; Nomoto et al. 2006). Through shocks formed by the ejecta colliding with 
the matter surrounding a SN, a fraction of this kinetic energy can be converted into X-ray 
emission: ~ 1% would be enough to explain the observation of SCP 06F6. 

The slow brightness rise, large photosphere size, the Swan spectrum, and strong X-ray 
emission set a general requirement for a supernova model: a dense C-rich environment around 
the progenitor. Let us estimate the density and mass in this envelope that are required to 
form a photosphere at Rph-eg ~ 60 AU . The optical depth in the envelope is: 

/■oo 

'T — Kpdr . (7) 

-Rph-eg 

For the photosphere t = |, and we will again use k ~ 0.1 cm^ g~^. Assuming a 1/i?^ density 
profile for the wind, the density of the envelope at the photospheric radius i?ph-eg becomes: 

p - T.4 . 10-- (^j^) " g^) " g c.-, (8, 

and the total mass of the envelope within radius it!ph-eg is: 

M„(fl,._.,) . 3 X 10- (sy^)"" (^)^4. (9) 

We can then estimate the wind mass loss rate required to get a mass of ~ O.OSM© within 
the 60 AU radius with a constant wind velocity, fwind: 

^ - (oT^) " (sfit) (i(&) ^« ■ d") 
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If the dense wind is concentrated towards a preferred plane and we see the object roughly 
edge-on, as suggested by the existence of both a cool photosphere and strong X-rays (see 
section E]), the involved mass can be somewhat lower. An outflow limited to ±30° from the 
plane would lower the required envelope mass and mass loss rate by a factor of 2. 

4-1.1. The SN Iln connection 

It should be noted that there exists a subclass of SNe interpreted as stellar explosions 
interacting with a dense circumstellar environment, the so-called type Iln SNe (Schlegel 
1990). Their spectra and lightcurves are strongly influenced by the circumstellar material 
(CSM) and some of them are strong radio and X-ray sources. It is noteworthy that of the two 
SNe that Gaensicke et al. (2008a) bring as roughly matching the lightcurve of SCP 06F6, 
SN 1994Y is type Iln (Ho et al. 2001) and SN 2006gy can be classified as a hybrid type 
between la and Iln (Ofek et al. 2007). 

SNe Iln are usually considered to be core collapse events, but they are a very hetero- 
geneous group (e.g., Filippenko 1997; Hoffman 2007) and it was suggested that they can 
also originate from thermonuclear SNe (Hamuy et al. 2003; Kotak et al. 2004; Aldering 
et al. 2006). The first case suggested to be a thermonuclear SN Iln was SN 2002ic, for 
which Hamuy et al. (2003) proposed a scenario involving an AGB companion in a SN la-like 
system. This SN became a prototype for the hybrid la/IIn type, also nicknamed Ha (Deng 
et al. 2004; Aldering et al. 2006) 

SNe Iln are characterized by narrow or moderately broad Balmer emissions superim- 
posed on the typical broad SN lines. For example, a few days after maximum of SN 2002ic 
mainly the narrow component of H^ line was visible (Hamuy et al. 2003). In this case, 
brightness in the I band was almost constant during 60d after maximum, but the lightcurve 
shapes among this group are very diverse. Another type Iln, SN 1994W, exhibited a 120d 
plateau after which the flux dropped suddenly, by 3.5 mag in the V band within 12 days and 
a similar amount in the R band (Chugai et al 2004). 

Typical estimated pre-SN mass loss rates of SN Iln are in the range 10^^ — 10^^ Mq yr^^ 
(van Dyk et al. 1996), while some of them are reported to reach even > 10^^ Mq yr^^ for 
a short time before the explosion (Chugai & Danziger 2003). The winds of the suggested 
la/IIn keep to the extreme side: 10~^ — 10~^ Mq yr~^ (Hamuy et al. 2003; Trundle et al. 
2008). For the controversial case 2006gy (Iln or la/IIn) the immediate pre-SN mass loss rate 
has been estimated at ~ 10~^ Mq yr~^ (Ofek et al. 2007) or even ~ 1 Mq yr~^ (Smith et 
al. 2008). On average, SN Iln are almost as bright in the optical as type la (Richardson 
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et al. 2002) and some recent examples are among the most energetic SNe ever observed - 
most of these have been proposed to be la/IIn (e.g., Smith et al. 2008). It has recently been 
argued that some core collapse SNe may explode from Luminous Blue Variables (LBV) as 
progenitors (Kotak & Vink 2006), in particular for type Iln (Gal- Yam et al. 2007). This 
led to the development of a modified SN-CSM interaction scenario for the brightest type Iln 
SNe, involving a very dense, opaque envelope close to the SN progenitor (Smith & McCray 
2007; Smith et al. 2008), formed by the intense LBV wind of the progenitor, most likely 
by a giant eruption of the rj Car type (as opposed to a companion wind in the SN la-like 
scenario of Hamuy et al. 2003). 

SCP 06F6, if located at z = 0.143, was optically bright, but not overluminous, Mj ~ —18 
(Gaensicke et al. 2008). Its X-ray brightness, on the other hand, was very high. A dense, 
optically thick C-rich circumstellar environment is required to explain its characteristics and 
there might be a hint of a narrow H„ line in the VLT and KECK spectra (notice a small 
spike at ~ 6560A in the assumed rest-frame spectrum of SCP 06F6 in Fig. [T]). Therefore 
SCP 06F6 could be a case similar to type Iln SNe, but with a C-rich envelope - while for 
typical Iln or even for la/IIn SNe the CSM is H-rich. 

Whatever the origin of the envelope, it has been shown by Smith & McCray (2007) 
that when it is opaque, photon diffusion of the thermal energy deposited in this envelope by 
shocked ejecta is able to reproduce the rounded and broad maximum light of SN 2006gy. This 
could be even more true for SCP 06F6, since the late rapid decline of this photon diffusion 
model was actually too fast for SN 2006gy - additional energy sources had to be included for 
the late light, which for SCP 06F6 would not be required. In the Smith & McCray (2007) 
scenario the faint X-ray emission of SN 2006gy is attributed to a later interaction of the 
shock with a less dense CSM, after the shock escapes the opaque envelope. The extremely 
strong but transient X-rays from SCP 06F6 differ from the SN 2006gy case and probably 
signify a different structuring of the CSM, e.g., a bipolar configuration. 

Other suggestions for the origin of SNe la/IIn include pair instability SNe (Smith et al. 
2007) and Quark Novae (Leahy & Ouyed 2008). These could also be tried for SCP 06F6, 
but exploring the more exotic possibilies is beyond the scope of the present paper. 

4- 1-2. Core collapse SN 

A core collapse SN is the explanation favored for SCP 06F6 by Gaensicke et al. (2008a). 
This is based on the spectral fit to a carbon star spectrum that gives z = 0.143, but the 
spectrum itself is not explained, instead they assume it to be a new class of SNe. A basic 
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problem with this scenario is that it requires a young massive progenitor, of a main sequence 
mass of > 8M0, while there is no bright star forming galaxy in the vicinity of SCP 06F6. 
There is a possible host galaxy at 1.5" from the transient position (Barbary et al. 2008) and 
at z=0.143, 1.5" is 3.7kpc. But if located at that redshift, this galaxy would be very faint, 
Mz 13.2 (Gaensicke et al. 2008a). 

Timescales of core collapse SN lightcurves are similar to that observed in the transient, 
but typically their declining parts are somewhat longer. There are many core collapse SNe 
which declined ~ 2 mag in the I band within lOOd. It is harder to find examples of ~ 3 mag 
decline within lOOd. Except for some SNe Iln, as explained above. 

The dense C-rich circumstellar envelope required by a SN model could be formed by a 
wind from the SN progenitor star itself or from its hypothetical companion. A companion 
could also naturally explain bipolar geometry of the matter. Winds from the canonical 
core collapse progenitors, Wolf-Rayet (WR) stars, can be dominated by carbon (the WC 
and WO types), but despite that C2 Swan bands were not observed in any SN. WR winds 
are characterized by fwind ^ 1000 km s~^ and M^ind ^ 4 x 10~^ M© yr~^ (Crowther 2007). 
Substituting these values in Eq. (fTOl) indicates that the density of a WR wind would be over 
3 orders of magnitude lower than required to form a photosphere at 60 AU. 

LBVs undergo giant eruptions during which their mass loss rate increases to IQ-^ Mq yr- 
or even ~ 1 M© yr~^ (e.g.. Smith Sz Owocki 2006). The high mass envelope could be pro- 
vided by such an eruption shortly preceding the SN event, as suggested by Smith et al. 
(2007) for the proposed LBV progenitors of SNe. However, for SCP 06F6 a problem arises 
that LBVs not only are not C-rich, but they are carbon depleted (e.g., Garcia-Segura et al. 
1996). Having both a WR and an erupting LBV in the system might provide the needed total 
density of matter and right abundance of carbon, but it would introduce an issue of wind 
interaction - it seems that the two winds would have to be mixed thoroughly to provide a 
carbon-dominated spectrum from the photosphere, lasting at least 36 days. Also, a chance 
to encounter a system with appropriate timing is small (Smith 2008 estimates that it may 
happen in 0.05% of all WR binaries). 

Finally, after the core collapse explosion, a few solar masses of SN ejecta moving with 
some 10"^ km s~^ would sweep up the few times 10~^Mq of the wind envelope contained up 
to the 60 AU radius in ~ 10 days, even if the latter was concentrated around the equato- 
rial/orbital plane. Therefore a preexisting steady wind cannot explain the relatively stable 
C-rich photosphere observed in SCP 06F6 at ~ lOOd from the start of the event. Only a well 
timed giant LBV eruption could supply enough matter to withstand the subsequent SN. 
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4.1.3. Thermonuclear SN 

Another possibility is that SCP 06F6 was a SN la-hke event occuring within an envelope 
of carbon-rich matter, a thermonuclear and C-rich version of SN la/IIn. The SN progenitor 
WD needs to accrete mass from a companion to exceed the Chandrasekhar mass. In this 
case the C-rich envelope would be explained by a carbon rich star being the mass donor to 
the accreting WD. The high mass loss rate demanded by Eq. (fTOl) suggests that the donor 
should be a cool AGB carbon star. For an AGB star the wind velocity is fwind ~ 10 km s~^, 
so the mass loss rate requirement from Eq. fllOl) would be ~ IO^'^Mq. This is the order 
of magnitude of the highest M„ind values determined observationally for AGB stars which 
are a few xlO^^M© (e.g., Justtanont et al. 2006). In cool carbon stars, the bulk of the 
envelope (and wind) is still hydrogen, but C2 bands are a prominent feature in their spectra 
up to ~ 5000 K (e.g., Barnbaum et al. 1996). If the photospheric temperature of SCP 06F6 
was significantly above that, the presence of Swan bands would imply carbon abundance 
somewhat higher than that typical for AGB carbon stars. 

As the explosion originates from an accreting WD, the ejection would likely form a 
bipolar structure from the very start. Also, the ejecta mass would be lower compared to the 
core collapse case. Both would make it easier for the wind equatorial belt to survive the SN 
blast and to feature a C-rich photosphere at 60 AU and ~ lOOd after the explosion. 

One can conjecture an even more exotic configuration, where the WD exceeds Chan- 
drasekhar's mass during a common envelope (CE) event with its AGB companion. The 
total envelope mass within 60 AU could then be comparable to the ejecta mass, making it 
even easier for the waist/disk structure in the envelope to withstand the SN blast and form 
a relatively stable photosphere of ~ 6000 K. A similar explanation was proposed for the 
prototypical la/IIn SN 2002ic by Livio & Riess (2003) who argued that it would favor the 
double-degenerate SN scenario. Ofek et al. (2007) suggest a CE connection to deal with the 
immense pre-SN mass loss rate of SN 2006gy. In any case, a binary companion in a natural 
way provides a bipolar geometry allowing for simultaneous cool photosphere of the donor 
wind and X-rays from the shocked bipolar ejecta (A different source of X-ray form a SN la 
was studied by Shigeyama et al. 1993). 

Normal SNe la have the IR lightcurve decline timescales similar to SCP 06F6 - in 
particular they can decline by ~ 3 mag in the IR bands within lOOd of the IR maximum, as 
did SCP 06F6. The initial brightness rise of the SNe is a few times faster, but in the scenario 
considered for the transient, it is the role of the dense envelope to slow down the initial 
lightcurve evolution. The maximum light peak would be smeared out by the interaction with 
the envelope, introducing the observed symmetry, the timescale (~ 200d), and the amount 
of brightness decline. Despite longer duration of the peak, additional energy deposited in the 
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opaque envelope by the trapped 7 photons and by the embedded coUision shocks can supply 
comparable (or even larger) amplitude of the optical peak as in a non-interacting case. The 
change from the photospheric to nebular phase of a SN will happen later than normal, and 
from that point on the lightcurve would fall down onto the radioactively powered exponential 
decline. The decline rate could resemble that of core collapse SNc because of the larger than 
normal for SN la total mass in the envelope, but the levels of the radioactive decay powered 
luminosity would be smaller because of the small mass of Ni ejected. 

The SN la-like scenario does not require a high mass progenitor, hence it can occur in a 
faint galaxy with no current star formation. The initial (main sequence) masses to produce a 
carbon star on the AGB lie in the range of ca. 1-4Mq (e.g. Lebzelter & Hron 2003; Stancliffe 
et al. 2005), which translates into main sequence lifetimes of a few xlO^ to 10^° yrs. 



4.2. IMBH-WD Tidal Disruptions 

A suggestion that SCP 06F6 might be a star disrupted by a BH was raised by Gaensicke 
et al. (2008; also Glcnnys R. Farrar, private communication 2008). The observed timescale 
of the event, the carbon-dominated spectrum, and the lack of a bright host galaxy lead us 
to consider a WD disruption by an IMBH, as recently modeled numerically by Rosswog 
et al. (2008b). We bring here a brief description of the tidal destruction of a WD by an 
IMBH. Future studies aiming at finding the light curve of such an event will have to take 
into account many other processes not considered here (e.g., Frolov et al. 1994; Wiggins & 
Lai 2000; Wilson, J. R. & Mathews, G. J., 2004; Brassart & Luminet 2008; Guillochon et 
al. 2008). 

A crucial parameter in star-BH collisions is the penetration factor, (5 — Ri/Rp, describ- 
ing how close the two objects get, i.e., how the pericenter of the orbit, Rp, compares to 
the tidal radius at which the tidal acceleration overcomes stellar (here - WD) self-gravity, 
= {Mbr/Mwd)^^^ RwD- Using typical CO WD mass of MwD = O.QMq and radius of 
-RwD = 9 X 10^ cm, and an IMBH mass of Mbh — lO^M©, the tidal radius is: 



Mbh \ ( ^WD \ / RwD 



Tidal disruption occurs for /9 > 1. But for a fixed /5, a WD will only be tidally disrupted by 
a BH within a certain mass range (Luminet & Pichon 1989a,b; Rosswog et al. 2008b). If the 
BH is too massive, the WD is swallowed before being tidally disrupted, with no significant 
emission. If the BH is too low-mass, it enters the WD instead of disrupting it. For too high 
values of the penetration factor there is no chance for a tidal disruption event at all, because 
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the strongly relativistic regime (the star enters the BH/the BH enters the star) would come 
into play at any mass. For intermediate masses and moderate values of /9, the interaction 
with the BH compresses the WD to a pancake-hke shape and tears it apart. Part of the 
material will be gravitationally bound to the IMBH and will form an accretion disk which 
can supply accretion luminosity close to the Eddington luminosity, mainly in soft X-rays, 
for about a year following the disruption (Rosswog ct al. 2008b), similar to the ~ 200 days 
duration of the SCP 06F6 event. Explosive thermonuclear ignition of the compressed WD 
material is also possible and Rosswog et al. (2008a) proposed this as a mechanism for a 
new type of underluminous SNe, later suggesting that SCP 06F6 can be the first discovered 
object of this kind (Rosswog et al. 2008b,c). But since SCP 06F6, if at 2; = 0.143, was clearly 
not underluminous compared to SNe, and accretion alone can supply enough luminosity, we 
do not consider a thermonuclear explosion as a necessary ingredient in this case. 

The optical observations do, however, require another element: enough matter around 
the object to form a relatively stable photosphere of ~ 6000 K at 60 AU. We postulate that 
this can be due to a strong wind from the accretion disk, colliding with the matter ejected 
at earlier phases of the disruption. Note that in this case the X-ray emission can come either 
from a shocked fast outflow material or from the inner parts of the accretion disk, in fact 
similar to what would occur in the WD-asteroid model. 

According to Luminet Sz Pichon (1989b), the maximum BH mass still allowing tidal 
disruption of a O.6M0 WD is AfBH,iim = 3 x 10^ M©. However, such a massive BH could only 
exert weak tidal effects without swallowing the WD. The optimal BH mass, allowing the 
strongest tidal effects, is MBH,opt = 1-53 x IO^Mq. For SCP 06F6, the minimum required BH 
mass can be estimated by comparing the BH Eddington luminosity to the SCP 06F6 lumi- 
nosity. The limit turns out to be Mbh ^ 2 x IO^Mq when considering only the optical peak of 
Lp ~ 1.3 X 10^2 ej-g s-^ with k ~ 0.1 cm^ g-^ as before, or Mrh > 2 x lO^M0(fi;/O.2 cm^ g-i) 
when including the X-ray emission (with k scaling appropriate for high temperatures). Ex- 
clusion of X-rays is justified if they originate in the post-shock gas of a fast outflow. With all 
these considerations, the BH mass range relevant for the IMBH- WD scenario for SCP 06F6 
is a few xlO^ to a few xlO^M©. 

Accounting for the Swan spectrum requires a considerable amount of carbon, either 
coming from the initial WD composition or produced during the interaction with the IMBH. 
The former is the case for the CO WD with /? = 1.5 from Rosswog et al.: after the disruption 
the material retains its initial composition (assumed in their simulation to be 50% oxygen and 
50% carbon). Even igniting CO WDs keep carbon-rich external layers. The latter could occur 
for a helium WD, when the tidal compression triggers a thermonuclear explosion. While for 
a strongly interacting He WD the final nucleosynthesis products are dominated by silicon 
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and iron (Rosswog 2008b), carbon can be the dominant end-product in a weaker interaction, 
when the compression phase is too short to fuse heavier elements (Luminet & Pichon 1989b). 
Note that such interrupted burning would produce an even more underluminous SN. 

In the standard disruption models, the mass accretion rate (and hence the accretion 
luminosity) falls off at late times as ~ t~^/^ or slower, where t = is set at the maximum 
of the interaction strength, much before the lightcurve reaches its maximum (Lodato et al. 
2008; Ramirez- Ruiz & Rosswog 2008). For SCP 06F6, there are very few measurements at 
late times, and their errorbars are quite large. However, it is clear that with t = at the 
beginning of the SCP 06F6 event, the late decline rate is faster than t~^/'^ both in the optical 
and in the X-rays. The X-ray signal dropped over two orders of magnitude within ~ lOOd. 
This again can be reconciled if a strong wind from the accretion disk depletes the disk in 
a few months and terminates the backflow onto the BH. Another possibility is that at late 
times a large fraction of the accreted energy goes into launching jets. 

The outer radius of the disk is -Rdisk — few x i?t (Rosswog et al. 2008b). This allows 
a simple estimate of the outflow terminal speed, influenced by the interaction with material 
residing around and still falling in: 

t^wind = e^^esc(i?disk) ^ 6 X 10^ {^'^^)) ^ ^^^^ 

where the scaling factor ^ encompasses our uncertainity of the disk size and of the wind slow- 
down. From Eq. ffTUj) with fwind = 6000 km s~^, the mass loss rate required for a photosphere 
at -Rph-eg = 60 AU is M„ind — 0.7 yr^^. Over the 200 days eruption, a total mass of a 
few O.IMq would be lost from the accretion disk and its vicinity. The outflow velocity and 
mass would resemble those of SN ejecta. 



5. SUMMARY 

We compared four scenarios for the SCP 06F6 transient event (Barbary et al. 2008), 
one galactic and three extragalactic (with a redshift oi z = 0.143 as suggested by Gaensicke 
et al. 2008a): 

1. A Galactic event of an asteroid colliding with a WD, 

2. A new class of core collapse SN, as suggested by Gaensicke et al. (2008a), 

3. Type la SN explosion, resembling the hybrid la/IIn type, 

4. Tidal disruption of a CO WD by an intermediate mass black hole (IMBH), as recently 
studied in detail by Rosswog et al. (2008b, c) and Lodato et al. (2008). 
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In principle, all four scenarios can lead to a transient event, and we expect that such 
events will eventually be discovered. However, as we presently attempt to explain SCP 06F6, 
we compare the models to the observations of this transient in Table 1 with the following 

comments (numbers as in the first column of the table): 

(1) Redshifted carbon star-spectrum does better than what we can find for a Galactic source. 
A Galactic source requires absorption Swan bands distorted by strong magnetic fields, possi- 
bly with addition of some emission Swan bands, which also puts its distance to be ~ 1.5 kpc. 
Despite its complexity, we cannot rule out such a Galactic source. The IMBH-WD disrup- 
tion model with no nuclear burning, i.e., a CO WD seems to do better in accounting for 
the spectrum. A SN embedded in C-rich matter is also fine, but the origin of carbon needs 
explanation. The la-hke SN scenario (SN la/IIn) seems more natural than core collapse, 
which would require a very unusual binary interaction. 

(2) For the Galactic source, total emitted energy is used to constrain the mass of the asteroid. 
For the extragalactic sources, the luminosity is within the range of SNe. For the IMBH-WD 
disruption model there is as yet no constraint from theory or observation on the luminosity. 

(3) As far as timescales, the SN models have the closest match, although the observed decline 
is faster than those of most core collapse SNe. Some SNe type Iln have equally fast decline 
in the IR bands and the SN la-like model (la/IIn) does good as well. The WD-asteroid 
model (where the event timescale is the Kelvin-Helmholtz timescale), were not observed in 
the past, but from theoretical point of view can explain the timescale. As for the IMBH-WD 
tidal disruption model, we must incorporate a new ingredient to existing calculations of the 
disruption process, which we suggest is a strong disk wind. 

(4) In the Galactic source the X-rays originate from close to the WD, as in dwarf novae. 
This can be the X-ray source at late times, when the space above and below the disk has 
been cleared from dense gas. Another possibility is the ejection of jets (or a collimated fast 
wind), which when coUide with the circum-system material become the X-ray source. The 
same two mechanisms can operate in the IMBH-WD tidal interaction model. It seems that 
the SNe models should also employ a non-spherical geometry that allows fast polar ejecta 
to be shocked and emit X-rays. This further disfavors a core collapse SN based on a single 
star model; we can safely conclude that if SCP 06F6 is a core collapse SN, its peculiarity 
comes from a strong binary interaction. We recall that the peculiarity in the SN la-like 
model comes from a carbon star with a high mass loss rate, that also transfers mass to the 
pre-exploding WD. Namely, in the SNe models of SCP 06F6 we attribute the peculiarity to 
a rare binary companion. 

(5) The non detection of a galaxy at the location of the event is more of a problem for core 
collapse SN models because the galaxy is fainter than known star-forming galaxies. Less of 
a problem to SN la-like, and even less of a problem for the IMBH-WD disruption model, 
as a galaxy with no recent star formation can do with these two models. Moreover, if we 
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Table 1: Comparison of the models 
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assume that the IMBH is the BH at the nucleus of the galaxy, then a very faint dwarf galaxy 
is indeed expected. The non detection is not a problem for the WD-asteroid merger model, 
as it is used to constrain the progenitor distance to -D > 1.2 kpc (Barbary et al. 2008). 

We summarize the strong and weak points of the four models in the last two rows 
of Table 1. The most important thing to note is that IMBH- WD tidal disruption model 
has weak points if we use the existing numerical simulations for this process as they are. 
However, we can solve these problems if we add a very strong disk outflow at velocities of 
several x 10^ km s~^. (i) Photosphere: With a mass of 0.1 — O.SM© being lost at 3000 — 
10^ km s~^, we can account for the photosphere (deflned at r = 2/3) at tens of AU, and 
in particular of 60 AU at maximum light, (ii) Lightcurve: With the rapid removal of mass 
from the disk, and the interaction of this wind with inflowing gas, we can account for the 
lightcurve declining rate that is much more rapid than t~^^^. This outflow is powered by 
accretion close to the Eddington limit, and possibly by extreme ampliflcation of magnetic 
fields in the accretion disk. 

For that, our favorite model is the IMBH- WD model. Second we rank the type la-like 
SN model, with a carbon star mass donor companion. Third we rank the WD-asteroid 
Galactic model. The core collapse model comes last, as it seems it must be very peculiar, 
and involve some kind of strong binary interaction. 

The predictions of the different models are straightforward, allowing for an observational 
test. The Galactic scenario implies that a faint WD, with a possible IR excess, resides at 
the event location. The extra-galactic scenarios predict a very faint galaxy in that direction. 
The type II SN model implies the presence of a recently active star forming region. 

We thank Glennys R. Farrar, Stephan Rosswog, Enrico Ramirez-Ruiz and Giuseppe 
Lodato for helpful discussions. This research was supported in part by a grant from the 
Israel Science Foundation, and from the the Asher Space Research Institute in the Technion. 
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